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Abstract 

We consider the regularized Ericksen model of an elastic bar on an elastic foundation on 
an interval with Dirichlet boundary conditions as a two-parameter bifurcation problem. 
We explore, using local bifurcation analysis and continuation methods, the structure of 
bifurcations from double zero eigenvalues. Our results provide evidence in support of 
Miiller's conjecture [18] concerning the symmetry of local minimizers of the associated 
energy functional and describe in detail the structure of the primary branch connections 
that occur in this problem. We give a reformulation of Miiller's conjecture and suggest 
two further conjectures based on the local analysis and numerical observations. We 
conclude by analysing a "loop" structure that characterizes {k, 3k) bifurcations. 

Keywords : microstructure, Lyapunov-Schmidt analysis, Ericksen bar model 
AMS subject classification: 34C14, 74N15,37M20 

1 Introduction 

In the late eighties J. M. Ball suggested that an interesting and important question in 
material science would be to understand the dynamical creation of microstructure [3]. A 
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model for creating microstructure would be a dynamical system with a Lyapunov functional 
that does not reach its infimum value on, say, the set of iyQ^'^(0, 1) functions, the infimum 
value being achieved instead by a gradient Young measure. Thus one might expect to 
obtain microstructure dynamically, hoping that as the Lyapunov functional decreases along 
trajectories, translates in time will form a minimizing sequence. 

One of the candidates for such a process proposed by Ball et al. in [1] is Ericksen's model 
of an elastic bar on an elastic foundation [6]. This is given by 

utt = {W\u,) + put,), - au (1.1) 

on the interval [0, 1] with the Dirichlet boundary conditions 

u(0,t) = n(l,t) = 0. (1.2) 

Here u := u{x, t) is the lateral displacement of the bar, (3 measures the strength of viscoelastic 
effects (the term (3uxxt provides the dissipation of energy mechanism) and a measures the 
strength of bonding of the bar to the substrate. W'{ux) is the (non-monotone) stress/strain 
relationship; in what follows we specifically take the double well potential 

W{z) = -y-lf. (1.3) 



It is easily checked that 
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\ul + 2W{ux) + dx (1.4) 



is a Lyapunov function for fll.ll) . 



Friesecke and McLeod [S] proved that fll.ll) admits an uncountable family of steady states 
that are energetically unstable but locally asymptotically stable. They also showed that 
initial data evolves, roughly, to a saw-tooth pattern with the same lap number, (i.e minimum 
number of non-overlapping intervals where the pattern is monotone) as the initial data. In 
other words, as Friesecke and McLeod put it in the title of their paper [7], dynamics is a 
mechanism preventing the formation of finer and finer microstructure. These results go some 
way to explain the earlier numerical results of Swart and Holmes |19j . 



MixUer [18] considered the regularized version of the Ericksen model, 

Utt = {W'{ux) + (3utx - lUxxx)x - oiu (1.5) 
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on the interval [0, 1] with the double Dirichlet boundary conditions 



m(0, t) = t) = 0; Ma;x(0, t) = t) = 0. 



(1.6) 



The main thrust of Miiller's sophisticated analysis was to describe the global minimizer of 
the associated energy functional, 



^ Jo <- ^ J 

Before we continue, we need to define periodicity more precisely. Consider a stationary 
solution uo of (II. 5p . Take its odd extension to [1,2] and identify the points x = and 
X = 2. If the resulting function is D2fc-periodic on the circle for some A; G Z, we say that 
uq is periodic. Then Miiller's result is that the global minimizer is a periodic function with 
a precisely defined dependence of the period on a and 7. He also suggested the following 
conjecture: 

Miiller's Conjecture |18] : Local minimizers of E2 are periodic. 

Recently, Yip [26] has proved this conjecture for solutions of small energy where W has the 
form 



In this case many calculations of energy of equilibria can be done explicitly. This case, with 
more general boundary conditions, was also considered in [20l [23] . Nucleation and ripening 
in the Ericksen problem with the above form of free energy density is considered from a more 
thermodynamical point of view by Huo and 1. Miiller 

Finally, in a related paper [23], an extension of Ericksen's model to system of two elastic bars 
coupled by springs as a model for martensitic phase transitions is mainly studied numerically. 

The dynamics of the regularized Ericksen bar (11.51) was investigated in [16], where global 
existence of solutions, existence of a compact attractor and convergence to equilibria was 
proved. Furthermore, the case of a = was investigated in detail and an almost complete 
characterization of the structure of the attractor was given in that case. A. Novick-Cohen 
has observed that if a = 0, the set of stationary solutions of (II. 5p is precisely the same as 
for the Cahn-Hilliard equation, which was thoroughly investigated in [TDl [H]; more work 
exploiting this connection between (II. 5p and the Cahn-Hilliard equation is in preparation |12j . 
In particular, the stationary solutions of (II. 5p for the double Dirichlet boundary conditions 




(1.7) 



W{p) = {\p\ - If. 
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correspond to the Cahn-Hilliard equation with mass zero. As a consequence the bifurcation 
diagram of the stationary solutions of fll.51) with a = contains only supercritical pitchfork 
bifurcations from the trivial solutions; only the branches without internal zeros can be stable. 

Other studies of the dynamics of (11.51) include the work of Vainchtein and co-workers 
[2H [22| [25] , who considered, in particular, time- dependent Dirichlet boundary conditions 
(loading/unloading cycles) in order to study hysteresis effects. 

In this paper we would like to 

• present some evidence towards verifying Miiller's conjecture and reformulate it; 

• explain how the situation for a = for (11.51) can be reconciled with the result of 
Friesecke and McLeod for (II. ip alluded to above. 

In very recent related work, Healey and Miller [13], have considered the two-dimensional 
version of the problem with hard loading on the boundary, using methods of global bifur- 
cation theory and numerical continuation techniques, concentrating on primary bifurcating 
branches and characterizing their symmetry. 

We use methods of local bifurcation theory. We start by obtaining the primary and secondary 
bifurcation points and presenting the results of numerical continuation using AUTO [5]. In 
section [3] we apply directly the Lyapunov-Schmidt theory as detailed in [5] ; this suggests 
a mechanism for the restabilization of unstable solutions. The analysis has uncovered an 
interesting pattern of primary branch connections which we analyse in section HI To conclude 
we present two further conjectures, these are based on the local analysis backed up with the 
numerical observations. 

2 Preliminaries 

We start by reviewing the bifurcation structure of the problem. As shown in [T^, the 
eigenvalues Uk of the linearization of (II. 5p around the trivial solution u = Ut = satisfy 

Pk = ^ (^pn^k^ ± ^/p^^l^k^ - 4{j7i^k^ - n^k^ + a)) . (2.8) 
Hence we have the following lemma. 
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Lemma 2.1 The eigenvalues of the linearization of U.5\) around the trivial solution u = 

Ut = are generically simple, and pass through zero at points 

n^kf — a , , 

7. = (2.9) 

where the integers ki are ordered by their distance from the number 



TT 

For example, if a = 25, ki = 2, k2 = 3, ks = 1, kn = n for n > 3. In other words at a = 25, 
based on the eigenf unctions sin^kirx), the first bifurcating solution branch has one internal 
zero, the second has two internal zeros, the third has no internal zeros and the nth has n — 1 
internal zeros. 

As we will be working in the (a, I/7) plane, it is convenient to use fl2.9p to define 

Tk = {{a, 1/7) I a e [0, 7r^k% I/7 = 'K^k^/{Ti^k^ - a)]. (2.10) 

Thus the curves are the curves on which the linearization has in its kernel the eigenfunction 
sui{kTTx). Note that F^ has the vertical line a = /c^vr^ as an asymptote. 

We can also determine the double zero eigenvalue points. The curves F^ and F; will intersect 
at a point where 



The corresponding values ■yk,i can be found from 

_ TT^fc^ - ak,i 

We call these bifurcation points (/c, /) bifurcations. Note that one can concoct any (/c, /) 
bifurcation point, but never a bifurcation point of multiplicity higher than two. From these 
double zero eigenvalue points curves of secondary bifurcation points emanate and we denote 
these by F^^. 



2.1 Numerical Evidence 

We use AUTO [5] to investigate numerically the bifurcation diagram of equilibria of fll.5p . 
that is we look at a, 7) = where 

a, 7) = -fUxxxx + au- {ul - Ux)x, u{0) = u{l) = 0, Uxx{0) = Uxx{l) = 0. (2.11) 
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For fixed a we can compute the bifurcation diagram in I/7 and two examples are sliown 
in Fig. [T] for a = 7.5 and a = 33. The figure shows bifurcations from the trivial solutions 
occurring at 7^^ and the secondary bifurcations {k, i), plotting the £^ norm of {u, Ux, u^xi u^xx) 
(an approximation of the norm) of the solution as 7 varies. For a = 7.5, {a)-{d) are 
the branches of solutions with 0-3 internal zeros and sample solutions on these branches 
are shown in Fig. [21 Sample solutions from the branches (e)-(/i) that bifurcate from these 
solution branches are also shown in Fig. [2l For a = 33 we have labeled the number of 
internal zeros for the branches that bifurcate from the trivial solution. 

a = 33 



L2 Norm 



a = 7.5 




(c) (f) 


' 


(e) 




1 1 1 





(g) 



100. I 150. 

(a) 



200. 250. 

l/y 




Figure 1: Two bifurcation diagrams for fixed values of a = 7.5 and a = 33. Figures show 
the £^ norm of (n, Ux, Uxx, Uxxx) as 7 varies. In Fig. [2] sample solutions are shown from each 
of the branches for a = 7.5; the arrows give an indication where each the solution is taken 
from. {a)-{d) show solution branches with 0, 1, 2 and 3 internal zeros. For a = 33 labels 1-7 
show the number of internal zeros. The "loop" for a = 7.5 connects solutions with and 2 
internal zeros whereas for a = 33 it is between solutions with 1 and 5 internal zeros. 



We can exploit the fact that the bifurcations {k, t) can be identified as limit points to perform 
two parameter continuation. The result of these computations is shown in the (a, I/7) plane 
in figure [31 We clearly see the curve Fi tending to the correct theoretical value of the 
asymptote at a = vr^. From the bifurcation (1, 2) the curve tends to infinite as a 
approaches zero as does Fsi; whereas the curves ri2, Fis, appear to tend to infinite for 
some a > 0. 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 




Figure 2: Sample solutions from the bifurcation diagram with a = 7.5. (a)— (d) 
shows solutions on the branches bifurcating from the trivial solution with I/7 = 
123.87, 111.34, 123.78, 218.2. (e)— (h) shows solutions on the secondary bifurcation branches 
with 1 /7 = 54.35, 121.40, 278.43, 226.78. 

3 (/c, /c + l)-bifurcations and a mechanism for restabi- 
lization 

To present a plausible scenario for restabilization of unstable equilibria, we are interested in 
the structure of stationary solutions of (12.111) in a neighbourhood of a. {k,k + 1) bifurcation 
point, when (i$(0,a,7) has a double zero eigenvalue with eigenfunctions = sin(/c-7rx) and 
ffc+i = sm{{k + l)7ix). To examine this we apply the Lyapunov-Schmidt theory as described 
in [9]. 

Set 

X = {ue C\iO, 1)), I n(0) = u{l) = 0, u^M = ^x.x.(l) = 0.}, 
and let Y = C°((0, 1)). We let L denote the hnearization : 

L = c?$(0,afc,fc+i,7fc,fc+i)- 

Let us examine the symmetries of (12.111) . Define on X two operators, -Ri and R2 by Riu = —u 
and R2u{x) = u{l — x). It is easily seen that the group {I , Ri, R2, R1R2} is isomorphic to 
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Figure 3: Two parameter continuation of the bifurcation points. For clarity only the curves 
Ffc = 1, 2, 3 and the secondary curves F^^^, /c = 1, 2, 3, ^ = 1, 2, 3 are plotted. 

Z2 © Z2, that $ commutes with this group, i.e. 

Ri^{u, a, 7) = ^{RiU, a, 7), 

and that if k is even, 

R2Vk = Vk while R2Vk+i = -Vk+i- 
Hence the theory of [9l Chapter X] is applicable. 
In the Lyapunov-Schmidt framework (see [9], Chapter VII]), 

X = sp{vk, Vk+i} © M, 
where M is the orthogonal complement to ker L. Similarly, 

Y = N ® range L, 

where N is the orthogonal complement of the range of L. Since L is self-adjoint, = kerL, 
so we take {vk, Vk+i) to be a basis of both the kernel of L and of A^. 
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By the theory of bifurcations with Z2 © Z2 symmetry, the bifurcation equation will be of the 
form 



where 





= 0, 




9iix,y,fi) 




Ax^ + Bxy'^ H 


- ajjLX 


92{x,y,fi) 




Cx^y + Dy"^ - 


f buy 



From now on we fix a at a point ak^k+ii and do not any longer indicate the dependence on 
a. We take our distinguished parameter to be A = I/7 and let it vary through the critical 
value l/7fc,fc+i. Clearly, 

d'^9i ^ d"^ d'^9i ^ 1 d^gi 

dXdx dX d'-fdx 7^ d'-fdx ' 

which means that the case of positive gi^x, say, corresponds to the case of negative giXx, so 

we are in case (A) of [9l p. 430]. Then varying a will unfold the degenerate bifurcation. 

By [HI Appendix 3] (see also (1.14) in P Chapter VII]), and since (f^ = by oddness, we 
get, for example, for gi, 

d^gi 



dx^ 

d^9i 
dxdy"^ 

d'91 
dxd'y 



= {vk, d^^{0,X){vk, Vk+1, Vk+i)) 
{vk, d^^{0,X){vk)), 



with similar expressions holding for the partial derivatives of g2- 
Note that, for example. 

Odd 

■\t,=t2=t3=0'i'[Vk 

i=l 

d d d 



d^^{Vk,Vk,Vk) - ^^T^\t 
Oti dt2 Ot-i 



i=l 

3 



3 



-l9,[{Vk)xf{Vk)xx- 



i=l 



Hence 



Similarly, 



A = - 



Id'gi 



6 dx^ 



3/c*7r'' / sin*(A;7rx) cos^(/c7rx) = -k^TT^ 
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D = ^(fc + l)V. 
9 



For a and h we have 



a = {vk, {vk)xxxx) = ^^^"^^ and h = {vt+i, {vk+i)xxxx) = ^{k + 1) V. 



Now to compute B and C: 



d d d 



d ^{Vk,Vk+l,Vk+l) = ■^-^-^\ti=t2=ts=0^{tlVk + {t2 + t3)Vk+l) 

- ^ ^ ^ \ti=t2=t3=0 [{Vk)xtl + {t2 + t3){Vk+l)x]l 



Hence 



B 



1 Q-'gi 

2 dxdy"^ 
»i 



dti dt2 dt3 

= -6 [iVk)xx[iVk+l)xY + '2{Vk)xiVk+l)xx{Vk+l)x] ■ 

3 [ k^{k + ifn^ sm^iknx) cos'^iik + l)7ix) dx 
Jo 



+ - / k{k + I)^7r^sm{2k7rx)sm{2{k + l)7rx)dx. 
2 Jo 

Since the second integral is zero, we have that 

B = C = 3k^{k + 1)2 [ sm\k7Tx) cos\{k + 1)ttx) dx = -7r\k + ifk^ 

Jo 4 

Now we can reduce the bifurcation equation to normal form. First note that 

ei = sgn (A) = 1, €2 = sgn (a)sgn (-1/7^) = -1, 

es = sgn (D) = 1, 64 = sgn (6)sgn (-1/7^) = -1, 
so that indeed we are in case (A) of [9l p. 430]. 

By Proposition 2.3 of [9l p. 424] the bifurcation diagram is determined by the modal pa- 
rameters 



m 



Da 



B and n 



Ah 



C. 



Now, 



and 



k^ 3(fc + l)4 4 
fc^ 8 13 



A;2 



n 



{k + iye = 2 



{k + lY?,k^A^ ' {k^lf 
Hence for all k, m > 2, but n can be either smaller or larger than one. If m > 1, n > 1, we 
are in region (1) of [H p. 433]; if m > 1, n < 1, we are in region (2). 
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This picture of local bifurcations from a double zero eigenvalue readily leads to a number of 
interesting conclusions. 

Let us first take the situation where both m and n are larger than one. The first such 
bifurcation is the (3,4) one. So assume that A; > 3. 

Consider the curves and F^+i in the {a, I/7) plane. For a < a^^k+i, Tfc lies below F^+i, and 
at the bifurcation point {ak,k+iA/lk,k+i) the situation reverses. The theory of Golubitsky 
and Schaeffer tells us that there are two branches of secondary bifurcations, which we call 
Ffc+i,fc and Tk,k+i, from the {k + l)-st and the k-th primary curves, respectively, such that 
close to the double eigenvalue point Tk,k+i lies to the left of the vertical line a = 0^,^+1 and 
Ffc_|_i lies to the right, see Fig. H] (a). These two curves of bifurcations form a wedge, which 
we call Kk- If below the wedge, to the left of ak,k+i the unstable manifold of solutions on 
the k-th primary branch has dimension and that of the {k + l)-st branch has dimension 
1, then to the right of aki under the wedge these dimensions switch, while in the interior of 
the wedge both the k-th and the {k + l)-st primary branches are stable. 



(a) (b) 




Figure 4: Local bifurcation structure as predicted by the Lyapunov-Schmidt analysis, in (a) 
m,n > 1 and in (h) m > l,n < 1. The shaded region indicates the wedge. 

In the case of m > 1, n < 1, using [9l Lemma 2.5, p. 426] we have that that both Ffc+i^^ and 
Tk,k+i lie to the left of the vertical line a = ak,k+i, with Ffc^^+i lying above Ffc+i^^; as before, 
this defines a wedge Kk in which both the k-th and the {k + l)-st branches are stable. This 
case is illustrated in Fig. H] (b) and confirms the numerical observation in Fig. [31 
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4 (/c, 3A;)-bifurcations 

During the process of numerical continuation of branches of solutions, we discovered an 
interesting "loop" structure connecting each k-th and 3/c-th primary branch (see Fig. [1]). 
We give a simple treatment based on Fourier mode truncations. We will investigate the 
(1,3) bifurcation in detail. By a scaling theorem of Aston [2], the structure of all {k,3k) 
bifurcations is the same. 

As the required computations are pretty involved, we proceed as follows: (1) we use a and 
7 as our parameters; this has the effect of changing hyperbolae into straight lines; (2) we 
work on the interval [0,7r], which allows us to get rid of many multiples of vr; and (3) we 
use MAPLE [17] throughout, particularly its polynomial manipulation abilities, such as the 
computation of resultants and of Grobner bases, and polynomial factorization. 

Thus, we are looking at 

3ulu^^ - u^x - 'ju^xxx - au = 0, xG[0,7r], (4.12) 
with double Dirichlet boundary conditions. 

Since we are close to the (1, 3) bifurcation point, we use the ansatz u = ai sin(x) +03 sin(3x). 
Multiplying equation (14.121) by each of sm{kx), /c = 1, 3 in turn, integrating from to tt and 
simplifying, we get the algebraic system 

1 82 1 3 3 ^ ,^ 

= --aai - -a^a-s + -ai - -a^ - 1/2701 
z o z o 

1 9 243 o 81 3 o 

= --aas + -as - —a^ - y 7^3 - g^i 

Let us examine (I4.13P in more detail. We need to work out what the curves Fi and F3 look 
like, and where they intersect. Putting 03 (ai) to zero in the first (second) of the equations 
of IKm gives that 

ri = {7 = l-a}, & F3 = {7 = 1/9 - l/81a}. 

These curves are plotted in Fig. [51 Thus in the (a, 7) plane above Fi there are no non-trivial 
solutions of (14.131) . We see that the (1,3) bifurcation point is at (a, 7) = (9/10, 1/10) where 
the two curves intersect. 



27 

T 
27 



alas 



(4.13) 
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Now let us see whether (14.131) will give us a loop, and what this loop really means. To have 
a loop we must have a solution with ai = and 03 7^ of equations fl4.13p . 

In other words, we need to solve simultaneously for the equations 

1 11 27 2 

— a H 7 Oq = 0, 

2 2 2 ^ 4 ^ 

and 

1 9 81 243 q ^ 

-2"«3 + - yT^s - -^^3 = 0- 

These solutions for 03 are the values of 03 on the pure branch at which a (pitchfork) 
bifurcation occurs. There are two equations in one unknown, 03, so for solvability they 
define a relation between a and 7. The form of this relation can be found by taking the 
resultant of the above two equations with respect to 03 and is given by 

7 1 

7 = a H . (4.14) 

' 153 17 ^ ^ 

This straight line of course intersects the lines Fi and r3 at (9/10, 1/10) and as the slope 
is positive, such bifurcation points only exist for a < 9/10. For example, it can be checked 
numerically that if a = 8/10, the (secondary) bifurcation from the 03 branch is at 7 = 
.09542483660. 

We now turn our attention to the nature of the "loop". For a < 9/10 there is no pure ai 
branch. The mixed mode branch bifurcates at 7 = 1 — a; the pure a-^ branch bifurcates off 
F3. It has a secondary bifurcation, a pitchfork, in the ai direction, on the line (I4.14p . One 
of these branches then hits the primary mixed branch and they disappear in a saddle node 
bifurcation. In fact, it is possible to work out where the turning points are. They are given 
by the following relation connecting a and 7: 

Pi (a, 7) :=8748 7 + 8748 a + 5557 - 2994732 aj^ 



+6016437 7^ - 1 19556 07 + 406782 7^ - 7938 
+1054782 0^7^ - 168492 0^7 + 243972 0^7 
-1000188 07^ - 3156 + 288684 7^ - 2187 = 0. 



(4.15) 



Equation (14.151) is obtained as follows: first divide the first of the equations of (I4.13P by ai 
and then compute the purely lexicographic Grobner basis of these two equations using the 
ordering 03 > ai. This results in a basis with two elements for the ideal generated by the two 
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original equations. One of these basis elements is found to be a polynomial in ai only. We 
can then take the resultant of this polynomial in ai with its derivative to obtain an enormous 
polynomial in a and 7, eliminating ai. Some of the roots of this polynomial correspond by 
construction to turning points, so it is just a matter of factorizing the resultant and picking 
the right term to obtain Pi (a, 7). If in Pi(a,7) we set 7 = 1/10 — 71, a = 9/10 — ai, we 
obtain a homogeneous polynomial of degree 4, 

P2(ai,7i) := 5557a^ - 2994732ai7i3 + 1054782a2^2 _ i68492af7i + 60164377^ 

If we now set 71 = hai, Descartes' rule of signs tells us that we will have two real solutions hi 
and h2 of P2(ai, hai)/a\ = 0. Using Maple to compute them, we obtain the approximations 
for the curves of turning points. 



7 ^ 1/10 - 0.203171(9/10 - a) and 7 ^ 1/10 - 0.043484(9/10 - a). 
(a) (b) 




Figure 5: (a) Plot in the plane of the bifurcations (b) plot illustrating the loop in ai, 
as and I/7 space. 



In Fig. [5] we plot in (a) Fi, F2, the pitchfork bifurcations given by fl4.14p and the curves 
of turning points as determined by fl4.15p local to the bifurcation point. We also illustrate 
schematically the bifurcations in 01,03,1/7 space. 

Remark: this phenomenon can also be analysed in the framework of bifurcations with 
hidden symmetries [H [TH [13] . 
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5 Two new conjectures 



We showed in section |3] that there are wedges Kk in parameter space (0,1/7) where both 
the kth and {k + \)st branches are stable. Unfortunately though this analysis is only local 
in nature. Local to the bifurcation point our numerics of section 12.11 agrees with the local 
analysis we have presented. This lends some confidence to the numerics. If we combine the 
local analysis with our observations from the numerical continuation in Fig. [3l we make the 
following conjectures. 

Conjecture 1. The wedges Ki^ are non-empty for all 7 < '^k,k+i- Along the curves Ffc+i fc, 
as a 0, 1/7 00, while along ^+1, I/7 ^ 00 as a ^ S/I 7^ 0. 

This automatically means that the curves Tk,k+i and Fj+i^/ intersect ii I > k. The upshot of 
this conjecture is that the different wedges Kk, Ki intersect, and increasingly so as 7 — > 0, 
which for each creates regions of parameter values in which there are different stable 
equilibria. This is partially consistent with the Friesecke and McLeod result, at least for 
a < Oi. Furthermore we can give a reformulation of Miiller's conjecture that we stated 
in the introduction, namely 

Miiller's conjecture. All stable equilibria are created by the above mechanism. 

The two statements of this conjecture are equivalent as all the stabilized branches have 
symmetry for some k. 

To conclude we propose another conjecture. 

Conjecture 2. If a > /c^vr^, there are no equilibria with less than [k — 1) internal zeroes. 

This is again based on the local analysis combined with the steady state bifurcation picture 
given in section 12. H the hope is this may be easier to prove than the other two. Note that 
if this conjecture is true, then for a > /c^vr^, an initial condition without internal zeroes will 
have to evolve at least k — 1 interfaces. 
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